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We review the assumptions and domain of applicability of Landau’s Hydrodynamical 
Model. By considering two models of particle production, pair production from strong 
electric fields and particle production in the linear o model, we demonstrate that many of 
Landau’s ideas are verified in explicit field theory calculations. 

PACS numbers: 


I. INTRODUCTION 

In the early 1970’s Peter Carruthers jj] [||] (at the 
prodding of one of his graduate students Minh Duong 
Van) realized that Landau’s Hydrodynamical model || 
explained the single particle rapidity distribution dN/dp 
of pions produced in the inclusive reaction P-P —> 7r 
+ X extremely well, (see Fig. 1) He then asked me to 
think about the effect of adding thermal fluctuations and 
with the help of Graham Frye and Edmond Schonberg 
at Yeshiva University (and also some help from Mitchell 
Feigenbaum) we found the correct covariant method of 
adding thermal fluctuations to the hydrodynamic flow 
which is now called the Cooper-Frye-Schonberg formula 
J|], H. Based on these ideas we were able to fit both the 
transverse momentum distribution as well as the rapidity 
distribution extremely well (see Fig. 2). 



FIG. 1. Comparison of the experimental C.M. rapidity 
distribution of outgoing pions from Proton Proton collisions 
at the ISR at (a) pisr = 15.4 GeV/c and (b) pisn = 26.7 
GeV/c compared with the no adjustable parameter result of 
the Landau Model from [js] 

The hydrodynamical model, which was considered 


“heretical” when applied to proton-proton collisions,as 
opposed to Feynman scaling models, was resurrected by 
Bjorken (j7j in 1983 to describe relativistic heavy ion col¬ 
lisions. The use of the hydrodynamical model to study 
heavy-ion collisions is now a sophisticated mini-industry 
using 3-dimensional fluid codes and sophisticated equa¬ 
tions of state based on lattice QCD. The effects of reso¬ 
nance decays have also been recently included. A recent 
fit to single particle inclusive spectra for Pb-Pb collisions 
at 160 AGeV and S-S collisions at 200AGev § is shown 
in Fig. 3. 



FIG. 2.Comparison of both transverse momentum and ra¬ 
pidity distributions for an ISR experiment with the Cooper- 
Frye-Schonberg reinterpretation of the Landau Model, from 

I 

After the initial successes of the hydrodynamical 
model, Carruthers and Zachariasen j9J |l(J made the first 
transport approach to particle production based on the 
covariant Wigner-transport equation. This further work 
of Pete’s inspired me and my collaborators David Sharp 
p[ and Mitchell Feigenbaum jlj] to study particle pro¬ 
duction in A0 4 in a mean field approximation using the 
formalism espoused by Carruthers and Zachariasen. Our 
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efforts in the 70’s were hampered by two factors: first 
we were not sure what the appropriate field theory was, 
and secondly computers were not large enough or fast 
enough to perform accurate simulations of time evolu¬ 
tion problems. This covariant transport theory approach 
turned out not to be the most convenient way to study 
initial value problems and would later be replaced by a 
non-covariant approach where the time was singled out. 



y 

FIG. 3. Fit of recent Pb+Pb and S+S data of the NA49 
and NA35 collaborations by Schlei et. al. 

Starting in the mid-80’s two new approaches were 
taken to understand time evolution problems in field the¬ 
ory, both related in spirit to the covariant transport ap¬ 
proach of Carruthers and Zachariasen which utilized a 
mean field approximation to close the coupled equations 
for the distribution functions. The first method was to 
assume a Gaussian ansatz in a time dependent variational 
method in the Schrodinger picture, [|ij| (l|. The second 
method was to directly study the time evolution of the 
Green’s functions in the Heisenberg picture in a leading 
order in large -N approximation [|l5[ [16| 0- 

Both these methods, which are mean-field approxima¬ 
tions, lead to a well posed initial value problem for the 
time evolution of a field theory. In the mean-field ap¬ 
proximation, we discovered that one has to solve at least 
10,000 equations for the Fourier modes of the quantum 
field theory in order to be in the continuum limit (i.e. 
for the coupling constant to run according to the contin¬ 
uum renormalization group). Thus, although much of the 
formalism was worked out by me and my collaborators 
in the mid and late 80’s, |l|] it was not until the ad¬ 
vent of parallel computation that numerical algorithms 
were fast enough to make these calculations practical. 
The first initial simulation attempts were presented in 
Santa Fe in 1990 at a workshop on intermittency that 


Pete asked me to help organize |lij. These first simula¬ 
tions took weeks of dedicated machine time. With the 
advent of the connection machine CM-5 at the ACL at 
Los Alamos simulations can now be done in a few hours 
finally making serious studies possible. 

In the past few years, we have been able to consider two 
aspects of Relativistic Heavy Ion Collisions. The first as¬ 
pect is connected with the production of the quark gluon 
plasma. The model we used was the popping of quark- 
antiquark pairs out of the vacuum due to the presence 
of Strong Chromoelectric fields. The mechanism we used 
was based on Schwinger’s calculation of pair production 
from strong Electric Fields [^0j. Using this model and as¬ 
suming boost invariant kinematics we were able to show 
that many of Landau’s assumptions were verified namely 
that from the energy flow alone one could determine the 
particle spectrum, and that the fluid rapidity spectra was 
the same as the particle rapidity spectra when one is in a 
scaling regime. Also the hydrodynamical prediction for 
the dependence of the entropy density as a function of 
the proper time was verified at long times. Numerical 
evidence for this will be displayed below. 

The second problem we considered was the dynamics 
of a non-equilibrium chiral phase transition. In this case 
because of the phase transition, one can obtain single par¬ 
ticle spectra which are different from a local equilibrium 
flow such as that given by Landau’s model. We found 
that when the evolution proceeds through the spinodal 
regime, where the effective mass becomes negative, low 
momentum modes grow exponentially for short periods 
of proper time. This leads to an enhancement of the low 
momentum spectrum over what would have been found 
in an equilibrium evolution. While considering this sec¬ 
ond problem we realized that the Cooper-Frye-Schonberg 
formula was valid even in a field theory evolution, pro¬ 
vided one interprets the field theory interpolating num¬ 
ber density as the single particle phase space distribution 
function of a classical transport theory |l7]] . 

It is safe to say that much of my career was stimulated 
by Pete encouraging me to understand multiparticle pro¬ 
duction in high energy collisions and by him freely shar¬ 
ing all of his intuition about this subject. It has taken 
20 years to show that our original thinking back in 1973 
was mostly correct!! 

Landau’s model || of multiparticle production was 
based on very few assumptions. He first made an assump¬ 
tion about the initial condition of the fluid. Namely, he 
assumed that after a high energy collision some substan¬ 
tial fraction (a: 1/2) of the kinetic energy in the center of 
mass frame was dumped into a Lorentz contracted disc 
with transverse size that of the smallest initial nuclei. He 
then assumed the flow of energy was describable by the 
relativistic hydrodynamics of a perfect fluid having an 
ultrarelativistic equation of state p = |.(This assump¬ 
tion was later modified by later workers as knowledge 
first from the bag model and then from lattice QCD be¬ 
came available). The motion of the fluid is described 
by a collective velocity field u^{x,t). Once the collec- 
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tive variables, energy density, pressure and collective ve¬ 
locity are postulated, along with the equation of state, 
the equations of motion follow from the conservation of 
energy-momentum. The same assumption of a perfect 
relativistic fluid is used in Cosmology. One solves the 
Einstein equations for the curvature assuming an energy 
momentum tensor that is a perfect fluid, (see @ ) 

The rest of the dynamics is embodied in the initial 
and final conditions. We stated already that the initial 
condition is to assume that the initial energy density dis¬ 
tribution is constant in a Lorentz contracted pancake. In 
our field theory calculations using strong fields, we will 
assume that all the energy density is in the initial semi- 
classical electric (chromoelectric) field, and equate the 
energy density E 2 (x,t) with Landau’s e(x,t). Thus we 
will not attempt to derive this initial condition but will 
retain this assumption. However various event generators 
do find similar energy densities to those obtained from a 
hydrodynamical scaling expansion. The final boundary 
condition is to state that when the energy density e re¬ 
duces to a critical value e c which is approximately one 
physical pion/ pion compton wavelength 3 in a cornov- 
ing frame, then there are no more interactions and one 
then calculates the spectrum at that “freeze out” surface. 
The original method of Landau to determine the parti¬ 
cle distribution was to identify the pion velocity with the 
collective velocity and assume that the number of parti¬ 
cles in a bin of particle rapidity was equal to the energy 
in that bin divided by the energy of a single pion having 
that rapidity. An alternative was to assume the number 
distribution was proportional to the entropy distribution. 
This ansatz was later modified by myself and my collab¬ 
orators by assuming there was a local thermal distribu¬ 
tion of pions in the comoving frame at temperature T c 
described by 


g{x , k ) = ^{expf/c^/Tj - 1 1 . 


The particle distribution was then given by: 


E d 3 N _ d 3 N 
d 3 k 7 T dk\dy 


g(x, tyk^dcTfi 


( 1 . 1 ) 


( 1 . 2 ) 


where is the surface defined by e = e c . We will ver¬ 
ify that both ideas can be justified by our field theory 
calculations. 

Let us now let us look at the hydrodynamics of a per¬ 
fect fluid. In the rest frame (comoving frame) of a perfect 
relativistic fluid the stress tensor has the form: 


T^ u = diagonal (e,p,p,p) (1.3) 


well as A-p and A-A collisions are a result of the hydro¬ 
dynamics being in a scaling regime for the longitudinal 
flow. (More exact 3-D numerical simulations with sophis¬ 
ticated equations of state have now been performed. The 
interested reader can see for example ||). 

That is for v = z/t (no size scale in the longitudinal 
dimension) the light cone variables r, p : 

2 = rsinh? 7 ;f = r cosh 77 (1-5) 

become the fluid proper time r = f(l — v 2 ) 1 / 2 and fluid 
rapidity: 

77 = l/ 21 n[(f — x)/(t + a;)] => l/ 21 n[(l — v)/(l + u)] = a 

( 1 . 6 ) 

Letting u° = cosh cc; rt 3 = sinha, we have when v = zt 
that p = a, the fluid rapidity. If one has an effective 
equation of state p = p(e) then one can formally define 
temperature and entropy as follows: 

f + P = 3-.;* = T<i»;ln» = /*/( e + p) (1.7) 

Then the equation: 

■u ,J d v T^ = 0 

becomes: 


<9"(s(t )u v ) = 0 (1.8) 

Which in 1 + 1 dimensions becomes 

ds/dr + s/t = 0 or st = constant (1.9) 

The assumption of Landau’s hydrodynamical model is 
that the two projectiles collide in the center of mass frame 
leaving a fixed fraction (about 1 / 2 ) of their energy in 
a a Lorentz contracted disc (with the leading particles 
going off). The initial energy density for the flow can 
then be related to the center of mass energy and the 
volume of the Lorentz contracted disk of energy. It is 
also assumed that the flow of energy is unaffected by the 
hadronization process and that the fluid rapidity can be 
identified in the out regime with particle rapidity. Thus 
after hadronization the number of pions found in a bin 
of fluid rapidity can be obtained from the energy in a 
bin of rapidity by dividing by the energy of a single pion 
having that rapidty. When the comoving energy density 
become of the order of e c we are in the out regime. This 
determines a surface defined by 

£ c{ T f) = m n /V 7: (1.10) 

On that surface of constant r, 


Boosting by the relativistic fluid velocity four vector 
u^(x, t) one has: 

T F = ( £ + p)«V-pf (1.4) 


dN 
dp 
da^ 
dN 
dp 


1 dE _ 1 r 

m^u, 0 dp cosh a/ dp 

Aj_(dz,—dt) = An a 2 Tf (cosh p, — sinlrz;) 

A± 

- -—\(e + p) coshacosh(?7 — a) — pcoshnl 

m v cosh a 


( 1 . 11 ) 


From a hydrodynamical point of view, flat rapidity 
distributions seen in multiparticle production in p-p as 
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where A± is the transverse size of the system at freezeout. 
If we are in the scaling regime where r] = a then 


dN 

dr] 


A i 


TO, 


<Tf). 


which is a flat distribution in fluid rapidity. At finite en¬ 
ergy, where scaling is not exact, only the central region 
is flat and instead one gets a distribution which is ap¬ 
proximately Gaussian in rapidity. Exact numerical sim¬ 
ulations (see for example (||) show that the isoenergy 
curves do indeed follow a constant r curve in the central 
region, so that the scaling result does apply for parti¬ 
cle production in the central rapidity region at high but 
finite center of mass energy. 

In Landau’s model one needed an extra assumption to 
identify the collective fluid rapidity a with particle ra¬ 
pidity y = 1/2 In [(En + p^)/(E v - p„)], where p „ is the 
longitudinal momentum of the pion. What results from 
our field theory simulations of both the production of a 
fermion-anti fermion pairs from strong Electric fields (the 
Schwinger mechanism [ 2(J ) as well as in the production of 
pions following a chiral phase tranisition in the a model, 
is that if we make the kinematical assumption that the 
quantum expectation values of measurables are solely a 
function of r we will obtain a flat rapidity distribution for 
the distribution of particles. We can prove, using a coor¬ 
dinate transformation, that the distribution of particles 
in fluid rapidity is exactly the same as the distribution 
of particles in particle rapidity. Furthermore we also w ill 
find that it is a good approximation to use eq. 1.11 to 
determine the spectra of particles. For the full single par¬ 
ticle distribution Ej^ we will find that the Cooper-Frye 
formula is valid with the identification of the interpolat¬ 
ing number density with the single particle distribution 
function of classical transport theory. 


II. PRODUCTION AND TIME EVOLUTION OF 
A QUARK-ANTIQUARK PLASMA 


which leads to a critical electric field of order 2n l c 

he 

The problem of pair production from a constant Elec¬ 
tric field (ignoring the back reaction) was studied by J. 
Schwinger in 1951 [^0| . The WKB argument is as fol¬ 
lows: One imagines an electron bound by a potential 
well of order | Vo | ~ 2 m and submitted to an additional 
electric potential eEx . The ionization probability is pro¬ 
portional to the WKB barrier penetration factor: 

r v o! e a 

exp[—2 / dx{2m(V 0 -\eE\x)} 1 ^ 2 }=ey:p(—-m 2 /\eE\) 

Jo ^ 

In his classic paper Schwinger was able to analytically 
solve for the effective Action in a constant background 
electric field and determine an exact pair production rate: 

00 /_p"+l 

w= [aE 2 /{ 27t 2 )]^---exp(— mrm 2 /\eE\). 

i ^ 

n=1 

By assuming this rate could be used when the Electric 
field was slowly varying in time, the first back reaction 
calculations were attempted using semi classical trans¬ 
port methods. Here we directly solve the field equations 
in the large-N approximation |lG|| . We assume for sim¬ 
plicity that the kinematics of ultrarelativistic high energy 
collisions results in a boost invariant dynamics in the lon¬ 
gitudinal (z) direction (here 2 corresponds to the axis of 
the initial collision) so that all expectation values are 
functions of the proper time r = \Jt 2 — z 1 .We introduce 
the light cone variables r and ??, which will be identified 
later with fluid proper time and rapidity . These coor¬ 
dinates are defined in terms of the ordinary lab-frame 
Minkowski time t and coordinate along the beam direc¬ 
tion 2 : by 


2 = Tsinh?7 , t = TCOshrj. (2.1) 

The Minkowski line element in these coordinates has the 
form 


Our model for the production of the quark-gluon 
plasma begins with the creation of a flux tube contain¬ 
ing a strong color electric held. If the energy density of 
the chromoelectric held gets high enough (see below) the 
quark-anti quark pairs can be popped out of the vacuum 
by the Schwinger mechanism |2C|j . For simplicity, here we 
discuss pair production (such as electron-positron pairs) 
from an abelian Electric Field and the subsequent quan¬ 
tum back-reaction on the Electric Field. The extension 
to quark anti-quark pairs produced from a chromoelec¬ 
tric held is straightforward. The physics of the problem 
can be understood for constant electric helds as a simple 
tunneling process. If the electric held can produce work 
of at least twice the rest mass of the pair in one compton 
wavelength, then the vacuum is unstable to tunnelling. 
This condition is: 

Tlj Q 

eE — > 2me 2 
me 


ds 2 = —dr 2 + dx 2 + dy 2 + T 2 dif . (2.2) 

Hence the metric tensor is given by 

9nv = diag(—1,1,1, t 2 ). (2.3) 

The QED action in curvilinear coordinates is: 

S= J d d+1 x (det U)[-l^V^ + |(Vt^) 7 ^ 

—im’F'F- (2.4) 

where 

= {dp + — ieA^ (2.5) 

Varying the action leads to the Heisenberg held equa¬ 
tion: 
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(7 M V M + m) T = 0, 


1 \ 7 3 

7 U ( d T + J + 7-L ’ <9± + - 1 — (d v - ieArj) + m 


and the Maxwell equation: E = E z (t) = —A v (t) 
= = 


( 2 . 6 ) 

7°7 3 Xs = A s x s 

( 2 . 12 ) 

= 0 , 

with \ s = 1 for s = 1,2 and X s = — 1 for s 
then get the mode equation: 

= 3,4, we 

(2.7) 

+UJ l- / ks (r) = 0 , 

(2.13) 

( 2 . 8 ) 

9 9 i 9 9 km CA 

4 + kl + m ; 7r„ = —-. 

T 

(2.14) 



The back-reaction equation in terms of the modes is 
= /[rfk](ki+m 2 )A 5 |/+| 2 , (2.15) 

S=1 J 

A typical proper time evolution of E and j is shown in 
fig. 4. Here an initial value of E = 4 was chosen. 


A. Spectrum of Particles 


To determine the number of particles produced one 
needs to introduce the adiabatic bases for the fields: 

^(®)= / Mk]^[6°(k;r)ii ks (r)e _l /“ k<ir 

** S 

+d( 0 )t (-k;r)u_ k 5 (r)e i / a;kdT ]e ik - x . (2.16) 

The operators b s ( k) and &i 0) (k; r) are related by a Bo- 
golyubov transformation: 


5(. 0 ) (k;r) = E a kr( r ) fe S ( k ) +/3 kr (r)4(-k) 

4° } (-k ;r) = 53^-(r)& 8 (k) + <(r)4(-k) (2.17) 


FIG. 4. Proper time evolution of E and j as a function of 
u = ln(^) for an initial E = 4. 


One finds that the interpolating phase space number den¬ 
sity for the number of particles (or antiparticles) present 
per unit phase space volume at time r is given by: 


We expand the fermion field in terms of Fourier modes 
at fixed proper time: r, 


*(*)= J [dk]J2lbs(k)^ s (r)e ikri e^ 

+4(-k)C ks (T)e-^e-* p ' x ] 


(2.9) 


The ^ then obey 


d_ _ 1 _ 

dr 2 r 


+ *7_l • kj_ + i^nrf + m 


n(k;r) = ^ (O in | 6 (, 0 )t (k; r)h(. 0 ) (k; r)| 0 in ) = E \(%r( T )f 

r= 1,2 s,r 

(2.18) 


Squaring the Dirac equation: 


V’ ks = 


d 1 \ . q 

— + — - * 7 ^ • k^ - *7 7 r^ + m 

clt 2 r 


This is an adiabatic invariant of the Hamiltonian dy¬ 
namics governing the time evolution of the one and two 
point functions, and is therefore the logical choice as the 
particle number operator. At r = tq it is equal to our 
± initial number operator. If at later times one reaches 

' ks ( T ) — 0 , ( 2 . 10 )^ ie 0U £ re gi me because of the decrease in energy density 
due to expansion it becomes the usual out state phase 
space number density. Although this does not happen 
for the above pair production in the Mean field approxi¬ 
mation, (because we have not allowed the electric field to 
dissipate due to the production of real photons), reaching 
(2 11 ) an ou t regime does happen in the a model if the energy 


> k s 

sR 
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density decreases as a result of an expansion into the 
vacuum. 

The phase space distribution of particles (or antipart- 
cles) in light cone variables is 


adiabatic phase space number density is one of two mo¬ 
menta canonical to the variables defined by the coordi¬ 
nate transformation to light cone variables. Namely the 
variables 


nk(r) = f{k v ,k±,T) 


d 6 N 

ti 2 dx^dk^dyd^ 


(2.19) 


A typical spectrum is shown in fig. 5 which 

shows the effect of the Pauli- exclusion principle. 
The raw results and also the results of averaging 
over typical experimental momentum bins are shown. 
This latter result compares well with a transport 
approach including Pauli-blocking effects (see [pT[). 


* 

a 



p - E -1 P 


FIG. 5. Comoving spectra of fermion pairs, before and after 
binning for an initial electric field E = 4 at r=400 


We now need to relate this quantity to the spectra of 
electrons and positrons produced by the strong electric 
field (the production of electrons and positrons from a 
strong electric field is our prototype model for the pro¬ 
duction of the quark gluon plasma from strong chromo¬ 
electric fields). We introduce the particle rapidity y and 
to_l = + m 2 defined by the particle 4-momentum 

in the center of mass coordinate system 

kfj, = (m± coshy, fcj_, m± sinhy) ( 2 . 20 ) 


r = (f 2 - z 2 ) 1 / 2 


77 = — In 
2 


t + z 
t — z 


have as their canonical momenta 


k T = Et/r — k z z/r k v = —Ez + tk z . 


(2.23) 


To show this we consider the metric ds 2 = dr 2 — r 2 dy 2 
and the free Lagrangian 


m dx^ dx u 
2 ds ds 

Then we obtain for example 

dr dt .dr dz. 


(2.24) 


Et — k,z 


= k^U,, 


(2.25) 


The interpolating phase-space density / of particles de¬ 
pends on k v , k_L, r, and is ^-independent. In order to 
obtain the physical particle rapidity and transverse mo¬ 
mentum distribution, we change variables from ( 77 , k^) to 
(z, y) at a fixed r where y is the particle rapidity defined 
by ( [2.20 ). We have 


E 


d 3 N 


d 3 N 


d 3 k Tidy dk 2 . 


= J tt dz dx\ J f(k v ,kj_,T) (2.26) 


where the Jacobian J is evaluated at a fixed proper time 
r 


dk v /dy dk^/dz 


dk v dy 

dy/dy dy/dz 


dy dz 

m± cosh (?7 — y) 

dk v 

1 


cosh 77 dz 


We also have 


(2.27) 


k T = m± cosh (77 — y); k v = —rm±_ sinh (?7 — y). (2.28) 


The boost that takes one from the center of mass coor¬ 
dinates to the comoving frame where the energy momen¬ 
tum tensor is diagonal is given by tanli 77 = v = z/t, so 
that one can define the “fluid” 4-velocity in the center of 
mass frame as 

= (cosliT?, 0 , 0 , sinh? 7 ) ( 2 - 21 ) 

We then find that the variable 



has the meaning of the energy of the particle in the co¬ 
moving frame. The momenta k v that enters into the 


Calling the integration over the transverse dimension the 
effective transverse size of the colliding ions A± we then 
obtain that: 

xdy dk\ =A± j dk r ' /(fcr? ’ k± ' T) “ 7rd?7 dk\ (2 ’ 29) 

This quantity is independent of y which is a consequence 
of the assumed boost invariance. Note that we have 
proven using the property of the Jacobean, that the dis¬ 
tribution of particles in partical rapidity is the same as 
the distribution of particles in fluid rapidity!! verifying 
that in the boost-invariant regime that Landau’s intu¬ 
ition was correct. 
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We now want to make contact with the Cooper- Frye 
Formalism. First we note that the interpolating number 
density depends on k v and k± only through the combi¬ 
nation: 


Wfc = Y + ~2 ~ (2.30) 

Thus f(k v ,k±) = f(k^u^) and so it depends on exactly 
the same variable as the comoving thermal distribution! 
We also have that a constant r surface (which is the freeze 
out surface of Landau) is parametrized as: 

da M = A±(dz, 0, 0, dt) = ^4j_dr?(cosh^, 0,0, sinh??) (2.31) 

We therefore find 



T 


k^da^ = Aj_m±T cosh(?y — y) = A±\dk v \ (2.32) 


FIG. 6. Time evolution of sr as a function of r 


Thus we can rewrite our expression for the held theory 
particle spectra as 


Here we have used 2.35 and the thermodynamic rela¬ 
tion: 


ndy dfc 2 = A± J dk vf( k ri’ k ±i T ) = J f u ^ T )k^ da ^ 

(2.33) 

where in the second integration we keep y and r fixed. 
Thus with the replacement of the thermal single particle 
distribution by the interpolating number operator, we get 
via the coordinate transformation to the center of mass 
frame the Cooper-Frye formula. 

Schwinger’s pair production mechanism leads to an En¬ 
ergy Momentum tensor which is diagonal in the(r, rj, x±) 
coordinate system which is thus a comoving one. In that 
system one has: 

T» u = diagonal {e(r),p||(r),pj.(r),pj.(r)} (2.34) 

We thus find in this approximation that there are two 
separate pressures, one in the longitudinal direction and 
one in the transverse direction which is quite different 
from the thermal equilibrium case. However only the 
longitudinal pressure enters into the “entropy” equation. 

Only the longitudinal pressure enters into the “en¬ 
tropy” equation 


£ + p || = Ts 


(2.35) 


d(er) 

dr 


+ P\\ 


Ejri 


d(sr) _ Ej eta 
dr T 

In the out regime we find as in the Landau Model 
st = constant 

as is seen in fig.6. 


de = Tds 

to calculate the entropy from the energy density and lon¬ 
gitudinal pressure. An alternative effective entropy can 
be determined from the diagonal part of the full density 
matrix in the adiabatic number basis. The energy density 
as a function of proper time is shown in fig.7. 



x 


FIG. 7. Time evolution of er as a function of r 


For our one-dimensional boost invariant flow we find 
that the energy in a bin of fluid rapidity is just: 


dE 

dij 


— J T°^da^ = A± t cosh t]e(t) 


(2.36) 


which is just the (1 + 1) dimensional hydrodynamical re¬ 
sult.Here however e is obtained by solving the field theory 
equation rather than using an ultrarelativistic equation 
of state. This result does not depend on any assumptions 
of thermalization. We can ask if we can directly calculate 
the particle rapidity distribution from the ansatz: 
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dN 1 dE A± 
c ??7 m cosh 77 dr] m 


(2.37) 


We see from fig. 8 . that this works well even in our case 
where we have ignored rescattering, so that one does not 
have an equilibrium equation of state. 



FIG. 8 . The ratio of the approximate rapidity distribution 
( to cosh y liy ) an( l exac t rapidity distribution as a function of 

T. 


In the field theory calculation the expectation value 
of the stress tensor must be renormalized since the elec¬ 
tric field undergoes charge renormalization. We can de¬ 
termine the two pressures and the energy density as a 
function of r. Explicitly we have in the fermion case. 

s(t) =< T tt >= tE s f [dk]R TT (k) + E 2 R /2 


where 

R T r(k) = 2 (pi + m 2 )(<?+|/ + | 2 - 5 0 ' l / _ | 2 ) - w 

— {p\ + m 2 )(n + eA ) 2 / ( 8 w 5 t 2 ) 

P\\( t )t 2 =< T vv >= tE s J[dk\A s 7rR vv (k) - ^E r t 2 (2.38) 


where 

R vrj {k)= 2|/ + | 2 — {2u)~ 1 {ui + A s 7t) - 1 — A s eA/8w 5 T 2 

— X s eE/8uj 5 — A s 7r/4o; 5 T 2 + 57 tA s (7t + eA) 2 /(16w 7 r 2 ) 

(2.39) 


and 

P±(t) = < Tyy >=< T xx A 

= (4r)' 1 ^2 f[ dk \{P±(P± + m 2 )~ 1 R TT ~ 2\Tip 2 L R riv } 

s J 

+E 2 R / 2. (2.40) 

Thus we are able to numerically determine the dynamical 
equation of state pi = Pi(e) as a function of r. A typical 
result is shown in fig. 9. 
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III. DYNAMICAL EVOLUTION OF A 
NON-EQUILIBRIUM CHIRAL PHASE 
TRANSITION 


Recently there has been a growing interest in the 
possibility of producing disoriented chiral condensates 
(DCC’s) in a high energy collision 122 123|.|24|] . This idea 


was first proposed to explain CENTAURO events in cos¬ 
mic r av e mperiments where there was a deficit of neutral 
pions ||25|l . It was proposed that a nonequilibrium chiral 
phase transition such as a quench might lead to regions 
of DCC ^4}. To see whether these ideas made sense 
we studied numerically pb|| the time evolution of pions 
produced following a heavy ion collision using the lin¬ 
ear sigma model, starting from the unbroken phase. The 
quenching (if present) in this model is due to the expan¬ 
sion of the initial Lorentz contracted energy density by 
free expansion into vacuum. Starting from an approxi¬ 
mate equilibrium configuration at an initial proper time 
r in the disordered phase we studied the transition to the 
ordered broken symmetry phase as the system expanded 
and cooled. We determined the proper time evolution 
of the effective pion mass, the order parameter < a > 
as well as the pion two point correlation function. We 
studied the phase space of initial conditions that lead 
to instabilities (exponentially growing long wave length 
modes) which can lead to disoriented chiral condensates. 

We showed that the expansion into the vacuum of 
the initial energy distribution led to rapid cooling. This 
caused the system, initially in quasi local thermal equi¬ 
librium to progress from the unbroken chiral symmetry 
phase to the broken symmetry phase vacuum. This ex¬ 
pansion is accompanied by the exponential growth of low 
momentum modes for short periods of proper time for a 
range of initial conditions. This exponential growth of 
long wave length modes is the mechanism for the produc¬ 
tion of disordered chiral condensates. Thus the produc¬ 
tion of DCC’s results in an enhancement of particle pro¬ 
duction in the low momentum domain. Whether such an 
instability occurs depends on the size of the initial fluctu- 





























ation from the initial thermal distribution. The relevant 
momenta for which this exponential growth occurs are 
the transverse momenta and the momenta k v = —Ez + tp 
conjugate to the fluid rapidity variable 77 = tanh _ 1 (z/t). 
We also found that the distribution of particles in these 
momenta had more length scales than found in local ther¬ 
mal equilibrium. When there is local thermal equilib¬ 
rium, the length scales are the mass of the pion and the 
temperature which is related to the changing energy den¬ 
sity, both of which depend on the proper time r. 

When we reexpress the number density in the comov¬ 
ing frame in terms of the physically measurable trans¬ 
verse distribution of particles in the collision center of 
mass frame, we find that there is a noticable distortion 
of the transverse spectrum, namely an enhancement of 
particles at low transverse momentum, when compared 
to a local equilibrium evolution. We will consider two 
cases, one in which there is exponential growth of low 
momentum modes due to the effective pion mass going 
negative during the expansion, and one where the initial 
fluctuations do not lead to this exponential growth. Both 
situations will be compared to a purely hydrodynamical 
boost invariant calculation based on local thermal equi¬ 
librium. In determining the actual spectra of secondaries, 
we find that the adiabatic number operator of our large¬ 
ly calculation replaces the relativistic phase space den¬ 
sity g(x,p) of classical transport theory in determining 
the distribution of pariticles in rapidity and transverse 
momentum. This makes it easy to compare our results 
with the hydrodynamical calculation in the boost invari¬ 
ant approximation which assumes the final pions are in 
local thermal equilibrium in the comoving frame. 

The model we use to discuss the chiral phase transition 
is the linear sigma model described by the Lagrangian: 


X = —\v 2 + A (a 2 + 7 T • 7 r) + XNGq(x, x). (3.5) 

If we assume boost invariant kinematics 1 0 which re¬ 
sult in flat rapidity distributions, then the expectation 
value of the energy density is only a function of the 
proper time. The natural coordinates for boost invari¬ 
ant (v = z/t) hydrodynamical flow are the fluid proper 
time r and the fluid rapidity 77 defined as 

Z t “h z 

To implement boost invariance we assume that mean (ex¬ 
pectation) values of the fields and x are functions of r 
only. We then get the equations: 

r~ 1 d T rd T $*(r) + x(r) $j(r) = HSu 

x{t) = A(—v 2 + $ 2 (t) + NG 0 (x , x■ r, r)), (3.6) 

To determine the Green’s function Gq{x,i)\t,t') we in¬ 
troduce the auxiliary quantum field <p(x, r) which obeys 
the sourceless equation: 

(r~ 1 d T rd T - r~ 2 d 2 - d\ + x(aO)</>(+r) = 0. (3.7) 

Go(x, y\ r, t') =< T{<j)(x,T) cfr(y,T')} > . 

We expand the quantum fields in an orthonormal basis: 

4>(ri,x±,T) = — J [d 3 k](exp(ffca:)/k(r) a k + h.c.) 

where kx = k v r] + kj_x± , [d 3 k] = dfc I? d 2 /cj_/(27r) 3 . The 
mode functions and x obey: 


L = • (9<E> — -A(<E> • $ — v 2 ) 2 + Ha. (3.1) 

2 4 

The mesons form an 0(4) vector *I> = (ct, 7 iy) This can be 
written in an alternative form by introducing the com¬ 
posite field: \ = A($ • $ — v 2 ). 

L 2 = ~(l>i{n + x)<t>i + ^ + ^XV 2 + Ha (3.2) 
The effective action to leading order in large N is given 

by @ 

T[$, x]=J d*x[L 2 ($, x, H) + l -Ntr lnGo 1 ] (3.3) 

Go X {x,y) = i[D + x(z)] S A {x-y) 

Varying the action we obtain: 

P + x{x)Vi =0 [□ + x(z)k = H, (3.4) 

where here and in what follows, TTi,a and x refer to ex¬ 
pectation values. Varying the action we obtain 


/k ~b (^2 + Ai( r ) + 4^2 )/ k — b- (3-8) 

X(r) = A(-v 2 + 4> 2 (t) + ^ n J [d 3 k]|/ k (r )| 2 (1 + 2 n k )). 

(3.9) 

We notice that when x goes negative, the low momen¬ 
tum modes with 


k 2 + 1/4 


+ k 2 ± < |xl 


grow exponentially. However these modes then feed back 
into the x equation and this exponential growth then 
gets damped. It is these growing modes that lead to the 
possiblity of growing domains of DCC’s as well as a mod¬ 
ification of the low momentum distribution of particles 
from a thermal one. The parameters of the model are 
fixed by physical data. The PCAC condition is 
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d n A ]i( x ) = Umlir l (x) = Hn l (x). 





In the vacuum state xo^o = 0 = H. so that <Jo = 

f n = 92.5 MeV. The vacuum gap equation is 


m 


2 

7T 


-Xv 2 + A fl + A N 



1 

2 v / fc 2 + ml 


This leads to the mass renormalized gap equation: 

X(r) -to 2 = -A/ 2 + A$ 2 (t) 

+ 7 W /[A]{|/ k( T )| 2 (It!*) - 2VF 1 +m , }. ( 3 . 10 ) 

A is chosen to fit low energy scattering data as discussed 



If we assume that the inititially (at tq = 1) the system 
is in local thermal equilibrium in a comoving frame we 
have 


n k 


1 

e^° E k - 1 


(3.11) 


where (3 0 = 1/T 0 and E% = +k\ + x(r 0 ). 

The initial value of x is determined by the equilib¬ 
rium gap equation for an initial temperature of 200-MeIA 
and is .7 /to - 2 and the initial value of o is just . The 
phase transition in this model occurs at a critical temper¬ 
ature of 160MeIG To get into the unstable domain, we 
then introduce fluctuations in the time derivative of the 
classical field. We varied the value of the initial proper 
time derivative of the sigma field expectation value and 
found that for tq = 1/to there is a narrow range of ini¬ 
tial values that lead to the growth of instabilities, namely 
.25 < \&\ < 1.3. 



T (fm) 


FIG. 10.Proper time evolution of the x field for two differ¬ 
ent initial values of a. 

Fig.10 displays the re sult s of the numerical simulation 
for the evolution of x ( |3.8|) (^T|). We display the aux¬ 
iliary field x i n units of frrT 2 , the classical fields <£> 
in units of /m -1 and the proper time in units of /to 
(1/to - 1 = 197MeF) for two simulations, one with an 
instability (<r| To = —1) and one without (<j| To = 0). 


We notice that for both initial conditions, the system 
eventually settles down to the broken symmetry vacuum 
result as a result of the expansion. The evolution of the 
quantities er and tt\ are displayed for various initial con¬ 
ditions in fig.ll. 





FIG. UProper time evolution of the <r field and n field for 
four different initial conditions with /, r = 92.5 MeV, l'25MeV. 

To determine the spectrum of particles we introduce 
the interpolating number density which is defined by ex¬ 
panding the fields in terms of mode functions f£ which 
are first order in an adiabatic expansion of the mode 
equation. 


p -iy k ( r ) 

fk = — ; dy k /dt = Uk, (3.12) 

V 2 Wfc 

where WkM = (^/f 2 + k'j_ +x( T )) 1 ^ 2 - This leads to the 
alternative expansion of the fields: 

(/>(r],x±,T) = J [d 3 k] (exp(ikx)(r ) a k (r) + h.c.) 

(3.13) 

The two sets of creation and annihilation operators are 
connected by a Bogoliubov transformation: 

afc(r) = a{k, r)a k + /3(k, r)aL fc - (3-14) 

a and (3 can be determined from the exact time evolving 
mode functions via: 
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and 


«M ^A) 

«*.r)=<(jJ^-^A). (3.15) 

In terms of the initial distribution of particles no(fc) and 
/3 we have: 

nk(r) = f(kri,k±,T ) =< a|.(r)afc(T) > 

= ?ro(fc) + |/3(fc,r)| 2 (l + 2?r 0 (fc)). (3.16) 


nk ( t ) is the adiabatic invariant interpolating phase space 
number density which becomes the actual particle num¬ 
ber density when interactions have ceased. When this 
happens the distribution of particles is 


f(k v ,k±,T ) 


d 6 N 

^dx^dk 2 ^ drjdkrj 


(3.17) 


We now need to relate this quantity to the physical spec¬ 
tra of particles measured in the lab. At late r our sys¬ 
tem relaxes to the vacuum and x becomes the square of 
the physical pion mass m 2 . As before, we introduce the 
outgoing pion particle rapidity y and m± = ^kj_ + m 2 
defined by the particle 4-momentum in the center of 
mass coordinate system. The boost that takes one from 
the center of mass coordinates to the comoving frame 
where the energy momentum tensor is diagonal is given 
by tanh ?7 = v = z/t, so that one can define the “fluid” 
4-velocity in the center of mass frame as 

= (cosh 77 ,0, 0, sinh?y) (3.18) 

The variable 


Wfc = Y + ^2 — (3.19) 


has the meaning of the energy of the particle in the co¬ 
moving frame. As discussed earlier the variables t, 77 have 
as their canonical momenta 


k T = Et/r — k z z/r k v = —Ez + tk z . (3.20) 

Changing variables from ( 77 ,/^) to ( z,y ) at a fixed r we 
have 


d 3 N d 3 N [ 2 T ,,, , , 

E #k = ^dkl = J ndz dx ± J /(fc - k± ’ T) 

= Aj dk v f(k vi k±, t) (3.21) 

This is again converted into the Cooper-Frye form by 
noting that a constant r surface is parametrized as: 

da M = A±(dz, 0, 0, dt) = Aj_d? 7 (cosh? 7 ,0,0, sinh 77 ) (3.22) 

Thus 


k^da^ = A±m±T cosh (?7 — y) = A±\dk v \ (3.23) 


-dy dk 2 = A± J dkr >f( k vi k -Li T ) = J f{k v ,k_ L ,T)k fl da ll 

(3.24) 

This reconfirms the idea that the interpolating phase 
space number density plays the role of a classical trans¬ 
port phase space density function, as was found in our 
calculation of pair production from strong electric fields 

We wish to compare our nonequilibrium calculation 
with the results of the hydrodynamical model in the same 
boost-invariant approximation. In the hydrodynamical 
model of heavy ion collisions 0 , the final spectra of pions 
is given by a combination of the fluid flow and a local 
thermal equilibrium distribution in the comoving frame. 

d 3 N d 3 N f 

<“ 5 > 

Here g(x, k) is the single particle relativistic phase space 
distribution function. When there is local thermal equi¬ 
librium of pions at a comoving temperature T c {r) one 
has 


g(x, k) = 3 w {exp[fc M u M /T c ] - 1} x . (3.26) 

In Figures 12 and 13 we compare the boost invari¬ 
ant hydrodynamical result for the transverse momentum 
distribution using critical temperatures of T c = 140, 200 
MeV to the two nonequilibrium cases represented in 
figure 10. Figure 12 pertains to the initial condition 
a\ To = —1 In this case there is a regime where the ef¬ 
fective mass becomes negative and we see a noticable 
enhancement of the low transverse momentum spectra. 
We have normalized both results to give the same total 
center of mass energy E cm . Figure 13 corresponds to the 
initial condition &\ To = 0. Here we notice that there is a 
little enhancement at low transverse momenta. 


da/dx=-1 initial conditions 



FIG. 12. Single particle transverse momentum distribu¬ 
tion for <7 = — 1 initial conditions compared to a local equi¬ 
librium Hydrodynamical calculation with boost invariance. 
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do/dx=0 initial conditions 
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FIG. 13. Single particle transverse momentum distribution 
for <7 = 0 initial conditions compared to a local equilibrium 
Hydrodynamical calculation with boost invariance. 

So we see that a non-equilibrium phase transition tak¬ 
ing place during a time evolving quark-gluon or hadronic 
plasma can lead to an enhancement of the low momen¬ 
tum transverse momentum distribution. In particular, 
if a Centauro type event is not accompanied by such an 
enhancement one would be suspicious of ascribing this 
event to the production of disoriented chiral condensates 
as a result of a rapid quench. For the case when we 
have massless goldstone pions in the a model (H = 0) 
P?! ,then x goes to zero at large times. In that case we 
find the amusing fact that the equation of state becomes 
p = e/3 at late times. This is true even though the final 
particle spectrum is far from thermal equilibrium. This 
equation of state is shown in Fig 14. 



) - ■ - 1 - ■ - 1 - ■ - 

0.0 50.0 100.0 150.0 


t 


FIG. 14. Equation of state | as a function of r for the 
massless a model where we start from a quench. 
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